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Abstract 

A model of scalar gluons and scalar quarks that successfully gives a 
Ins behavior in high energy qq scattering and which contains a three 
gluon vertex is used to derive transport equations for the quarks and 
gluons. Quasiparticle and semiclassical approximations are utilized. In 
particular, the collision integral is studied. At first the terms are orga- 
nized according to the number of gluon exchange interactions, and it is 
shown explicitly that the occurrence of all two body differential scat- 
tering cross sections in the Boltzmann equation requires an evaluation 
of all terms to this order. One obtains qg — > qg, qq — > gg, qq qq, 
and qq — >■ qq. An additional ordering according to the inverse num- 
ber of colors reveals that the quark loop diagrams are suppressed, and 
diagrams that lead to gluon production or scattering dominate. The 
self energy with three gluons exchanged is then examined. From this, 
it is evident that initial as well as final state interactions enter into a 
generalized collision integral. The generalized Boltzmann-like equation 
involving n gluon exchange for quark transport is given. 
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I. INTRODUCTION 



A central issue in describing heavy ion collisions is the question of particle frag- 
mentation, for which there exist several models fl]]. From a theoretical point of 
view, a description of such a collision usually is made within the framework of 
non-equilibrium theory in which the equation governing the Wigner function that 
describes the particle distribution f{X,p) semiclassically must somehow be given. 
From our knowledge of classical systems out of equilibrium and the attendant Boltz- 
mann equations, it is possible to make heuristic ansdtze for the equations governing 
the distribution functions of the relevant degrees of freedom. This line of thought 
has been taken up by several authors in particular from the partonic point of view 
0, and, based on this, parton cascade codes have been set up that are successful 
and necessary for the description of hard processes. 

On the other hand, a convincing formal justification of these heuristic ideas using 
Green function techniques does not yet exist, although steps have been taken in this 
direction in attempts to derive transport equations directly from the underlying 
theory of the strong interactions, quantum chromodynamics (QCD), see for example 
0-^. By now, the basic equations of motion for Green functions within a non- 
equilibrium approach for any theory are well known: they can simply be written 
down in a formal sense; one has a complicated set of matrix Dyson equations relating 
matrices of Green functions with matrices of self energies, for each particle species. 
The crux of the issue, after Wigner transformation, is the determination of these 
self energies for a given theory. The lowest order or mean field contributions to the 
self energy has been studied in several models of bosons or fermions or both, see for 
example [|TU],^. Higher order terms that contribute to the self energy and which 
can be interpreted as giving rise to collisional effects are less well studied. Recently, 
however, it was shown exactly within a simple relativistic model of quarks, the 
Nambu-Jona-Lasinio model [Q, that all graphs that contribute to second order in 
the interaction are required to generate two body elastic scattering diagrams that 
contribute to the differential cross-sections occurring in the Boltzmann equation 



In this current paper, we study a partonic model of QCD inspired by Polk- 
inghorne JI^ that contains scalar quarks and scalar gluons that are coupled to each 
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other. In addition, the gluons interact with each other via a triple gluon vertex, 
simulating the non-abelian nature of QCD. Such a 3-gluon vertex turns out to be 
sufficient to produce "pomeron-like" behavior i.e. the quark-quark elastic scattering 
amplitude can be calculated in a perturbative expansion in leading logarithms of the 
center of mass energy s, and such an expansion involves only gluonic ladder like dia- 
grams constructed between the two scattering quarks [^]. As it turns out, a 4-gluon 
vertex would be subleading in the expansion in In s and is therefore omitted in this 
model Lagrangian. We will take this as our starting point. It has the advantage over 
QCD proper that we can study the consequences of the non-abelian like terms, but 
we do not have to concern ourselves with the additional problem of gauge invariance 
|]6|,p!6[|. There is one additional difficulty, however, that the color structure of the 
theory becomes necessarily more complicated, but this can be handled relatively 
easily. 

Using the scalar partonic model, we examine the self energies that are required in 
the non-equilibrium formalism in and beyond the mean field approximation for use 
in constructing a collision term. Within the mean field approximation, the Hartree 
term is identically zero, due to color traces. The Fock term, which is often ignored 
in such calculation^, is discussed in some detail. The assumption that this term 
vanishes also, leads to enormous simplifications in considering higher order terms in 
the coupling g. We then examine the two gluon exchange diagrams, of which there 
are several different kinds. We classify them according to their topology as being of 
the "rainbow", "cloud", "ladder", "exchange" or "quark loop" kind. We find that 
all the graphs occurring in this order are required in order to produce the two body 
partonic differential scattering cross sections that would enter into a Boltzmann-like 
transport equation, in contrast to that suggested in 0]. This is an exact result. In 
our case, regarding for example the equations of motion for the quark distribution 
function, we find that the relevant cross sections that enter into the collision term at 
this level are either purely quark like in nature {qq qq and qq qq), or contain 
gluon components {qg qg and qq gg)^- Due to the plethora of different types 



^See however Refs |17| and [18| 



^ Other processes, such as qqg g and qqgg —>■ are suppressed by kinematical 
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of graphs that already occur at the two gluon exchange level, we seek an additional 
expansion parameter to simplify this and to enable us to go to higher orders in 
the coupling. In QCD itself, an expansion in the inverse number of colors is often 
useful, and we thus examine such an expansion within this model. An analysis of 
the color structure of each graph naturally suppresses diagrams containing quark 
loops and the exchange graphs thereof, as is also a feature of QCD. All remaining 
self energy graphs that contain gluonic exchange components (the ladder, rainbow 
and cloud graphs) turn out to be superficially of the same order in Nc, and all are 
required to generate all possible types of gluon production or quark-gluon scattering. 
A closer look, however, shows that the ratio of the ladder to the rainbow or cloud 
graphs goes as the ratio of the Casimir operators of the adjoint to the fundamental 
representations squared, which favors the ladder graphs numerically. The two body 
process is then briefly reanalysed along these lines. 

With this point of view, the next order in the coupling, i.e. the three gluon ex- 
change is investigated. Examining once again the equation of motion for the quark 
distribution function, we find that three parton production, i.e. qq —>■ ggg and 
(19 ^ Q99i naturally occurs. In addition to these production processes, however, 
the processes qq — > gg and qgg qg must also be accounted forQ. We thus see 
that initial as well as final state processes must be included in the matrix elements. 
With this information, the generalized Boltzmann equation that includes multipar- 
ton production {qq ng, but also processes of the form qqmg — > ng, with m and n 
integer) within this scalar quark and gluon model is given. 

It is to be noted that, unlike the situation that occurs in evaluating the elastic 



quark-quark scattering amplitude within the scalar parton model |T^, we cannot a 
priori make an assumption on the internal ordering of momenta within a self energy 
graph. This comes about in the elastic quark-quark scattering amplitude, because 
there are two external momenta that are fixed. In a self energy graph, however, 
and in the Boltzmann equation, only one external momentum can be fixed while 



considerations. 

^Again, such processes as qqggg — > and qqgg —>■ g are suppressed by kinematical 
considerations. 
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all the remaining momenta are integrated over. Thus the simphfications leading 
to the ladder approximation and Regge behavior in that calculation are not simply 
recovered in the kinetic theory. 

This paper is organized as follows. In Section 0, we give a brief introduction to 
the scalar partonic model, and set our notation for the nonequilibrium dynamics. 
In Section |I|, we describe the collisional contributions that arize from the mean 



field self energies, and the associated transport equations. In Section 0, we deal 
with the collision integral beyond the mean field approach, firstly giving the exact 
perturbative calculations with two gluon exchange, the color expansion, then three 
and n gluon exchange. We summarize and conclude in Section 0. 



II. MODEL OF SCALAR QUARKS AND GLUONS, AND TRANSPORT 

THEORY 

In this section, we introduce and discuss the scalar partonic model, and give the 
equations of motion for quark and gluon fields. We then give a brief discussion of 
transport theory applied to this scalar model, and indicate the self energies that we 
are required to evaluate. 



A. Scalar partonic model 



In the model underpinning our calculations, quarks and antiquarks are described 
by complex scalar fields (p, and gluons as the scalar field x coupled through the 
Lagrangian 

1 7n ^ 

£ = d^^^'%^i,l + -d^Xa^rd.x"^' - —Xa,rr'' 

-gm<P^^\Ty,{TnT<Pj,mXa,r - ^-^fatcfrstX^''x'''x'''. (2.1) 

The quark fields are regarded as massless, as one generally assumes for high 
energy processes, while the gluons are usually assigned a mass m a priori in order 
to avoid infra-red divergences. There is an interaction between quarks and gluons 
as well as a cubic self-interaction between gluons. Since in QCD the quartic inter- 
action between gluons leads to terms which are sub-leading in In s, such a quartic 
interaction among gluons is not included within this model. 
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One notes that both the quark and gluon fields carry two labels. Both of these 
refer to color groups. The fact that a direct product of two color groups is necessary 
can be seen on examining the three gluon vertex term. This term must be symmetric 
under the exchange of two gluons since they are bosons. In addition, one expects 
that the interaction vertex should be proportional to the (antisymmetric) structure 
constants of the color group. A single color group cannot meet these requirements, 
and the simplest combination which can is a product of two SU [Nc] groups. Thus 
the gluon field carries two color indices (a, r = 1...{N^ ~ !))• Since the quark field 
transforms in the fundamental representation of both of these SU{Nc) groups, it 
must carry two color indices as well = l...Nc). The matrices and are the 
generators and fate and frst are the structure constants for the two SU{N) groups 
respectively. Thus they satisfy 

[T^ T'] = ifabcT\ [T^ T^] = ifrstT'. (2.2) 



Note in Eq. ( p.l| ) that the fiavor index of the quark fields is suppressed. 



The equations of motion for the fields can be derived from the Euler-Lagrange 
equations. They are 

= -gm{T%{T^)l<P^'^^^'^Xa,r (2.3) 

for the (anti-)quarks and 

(□ + ^2)^a,. ^ _^^[0t«.'(r'^)5(r^)-0 . + r'^r'xi.,sXcA (2.4) 

for the gluons. 



B. Transport theory for interacting bosonic fields 

For the purpose of establishing our notation, we give the basic definitions and 
refer the reader to standard texts [0-0]. The quark Green functions in the Schwinger- 
Keldysh formalism [|l9l are defined as 

iS\x,y) = (T0^''(x)0ti.'^(y)) - (0^''(x)) ((j)^^^"'{y)) = tS—{x,y) 

iS\x,y) = (T(j)''\x)(f)^^'"'{y)) - (0*''(x)) ((^^^^"'{y)) = tS++{x,y) 

iS>{x,y) = (0*''(x)0tJ''»(y)^ - ^0*''(x)) (0^^'™(y)) = iS+-{x,y) 

iS<{x,y)= (0ti,'n(^)0M(a;)^ _ ^0i,«(a;)^^0ti,m(^)J) =is-^{x,y) (2.5) 
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and for the gluons as 

iG\x,y) = (Tx'^'^(a;)/'^(y)) - (x'^'^(x)) {x'^^y)) = iG" {x,y) 

iG\x,y) = (Tx"'^(x)x''1y)) - {X^^'^ix)) {x'^'iy)) = iG+^ix.y) 
iG>{x,y)= =tG+-{x,y) 
iG<{x,y)= = iG-^{x,y). (2.6) 

Here T and T are the usual time and anti-time ordering operators respectively. 
As given, the Green functions fall along the contour designated in Fig. 1. Our 
convention follows that of Ref. They satisfy a Dyson equation that introduces 
the matrix of self energies for either the quark or gluonic sectors, or S^. Using a 
generic notation, D_ = S_ot G and n = or as appropriate, one may write 

D{x, y) = D^x, y)- j d^zd^wD^x, w)n(w, z)D{z, y) 
= D^ix, y)- j d'^zd^wD{x, w)Il{w, z)D%z, y). 

(2.7) 

In a standard fashion, the transport and constraint equations can be derived, and 
this is summarized briefly in Appendix A. In terms of the center of mass variable 
X = X + y and the Fourier transform variable p conjugate to the relative distance 
u = X — y, or Wigner transform, the equations that one obtains for the off diagonal 
Green functions, corresponding to the transport and constraint equations are 

-2ipdxD-^{X,p) = I_ transport (2.8) 

and 

Qdx - 2p2 + 2M^^ D-+{X,p) = J+, constraint (2.9) 

respectively, where M is a generic parton mass, M = m for the gluons and M = 
for the quarks. is an abbreviation for the combined functions 

= JeoU + + /^?, (2.10) 

and /coll is the collision term, 

/,oU = U-+{X,p)AD+~{X,p) - U+-{X,p)AD~+{X,p) 

= Cr-/coU- (2.11) 
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and are terms containing retarded and advanced components respectively, 



= -U-+{X,p)AD^{X,p) ± D^{X,p)AU-+{X,p) (2.12) 

and 

= n^(X,p)AD-+(X,j9) T D-+{X,p)AU''{X,p). (2.13) 

In Eqs. (p.ll| ) to (|2.13|) , the operator A is given by 

A := exp 1^ (dx~^p - tptx)]. (2.14) 

In order to proceed further, we have to calculate the self energies. It is useful to 
introduce the quasiparticle approximation, in which a free scalar parton of mass M 
is assigned the Green functions 

zD-+(X,p) = ^{5{E, - p')fa{X,p) + 5{E, + p')f-a{X, -p)} (2.15) 
tD+-{X,p) = ^{6{E, - p')fa{X,p) + 6{E.^+p')MX, -p)} (2.16) 

tD-{x,p) = , + e{-p')D^-{x,p) + e(/)D-+(x,p) 

+ ^{6{E, - p')fa{X, p) + 5{E, + p°)/.(X, -p)} 

(2.17) 

zz}++(x,p) = ^.^ + 0( + e(/)D-+(x,p) 



p2 


-M2 + 






z 




p2 








— i 




p2 


-M2- 


ie 








p2 


-M2 - 


ie 



+ —{6iE, - p')faiX,p) + 5iEp+p")MX, -p)} 

(2.18) 

with Ep = p"^ + M^, and which are given in terms of the corresponding scalar quark 
and gluon distribution function, fa{X,p), and fa = 1 + /«, where a denotes the 
parton type a = q,g. 

Our task in this paper is to construct an equation for the distribution functions 



for quarks and gluons fa{X,p) from Eqs.(|2.8|) to ( p.l3|) , using the quasiparticle Green 
functions of the form given in Eqs.( |2.15D to ( p.l8|) . To do so, it is necessary to inte- 
grate the entire Eq.( |2.8| ) over an interval A± which contains ±Ep{X). Our particular 
focus will be placed on the collision terms that are derived from Eq.(^?TT]). In an- 
ticipation of this fact, we will require such an integration over the terms appearing 
in Eq. (|2.11| ), and define 
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T _ Tgain _ tIoss 
■^coll — <^coll ■^coU 

= Idpo^S- [dpol^t- (2.19) 
ja+ ja+ 

To lowest order in an expansion that sets A = 1, this integral can be easily performed. 
One has 

Jcoii=/ dpoU-+{X,p)D+-{X,p)- f dpoU+-{X,p)D-+{X,p) 
ja+ Ja+ 

= n-+(X,po = E„p) UX,p) + n+-(X,po = E,,p) fa{X,p), 

hip hip 

(2.20) 

i. e. the off-diagonal quasiparticle self energies are required to be calculated on shell. 

III. THE COLLISION INTEGRAL - MEAN FIELD SELF ENERGIES 

In this section, we evaluate the self energies to first order in the interaction 
strength and illustrate their role in the transport equation in the semi-classical 
limit. 

A. Hartree self energies 

For the scalar parton model, two generic kinds of Hartree graphs can be identified 
in the quark and gluon self energies. These are depicted in Fig. 2. For Hartree 
diagrams of any kind, off diagonal self energies are per definition zero and only 
diagonal elements can possibly be constructed, i.e. or S^"^. However, all such 
diagrams vanish identically in this model. The reason for this lies in the color factors: 
for the quark self energy graph in Fig. 2 that contains a quark loop, a single SU(3) 
color group leads to the associated color factor 

Fh,, = tli^m = 0, (3.1) 

since ta — A"/2, where A" are the Gell-Mann matrices. In the above expression, 
i denotes the external quark momentum and is therefore not to be summed over. 
For the quark self energy containing a gluon line, the color factor for a single SU(3) 
group is also vanishing. 
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FH,g - tiiT^i, - -itiJabb - 0. (3.2) 

Similar arguments apply to the gluon self energies. Thus, no mass renormalization 
occurs due to Hartree terms. 

A semiclassical expansion of the transport and constraint equations thus leads 
to free streaming described by 

-2ipdxD-+{X,p) = 0, (3.3) 

and the constraint equation is 

{p^-M^)D-+{X,p) = 0. (3.4) 

By integrating Eq.( p.3|) over a positive interval A+, the equations for the time evo- 
lution of the distribution function for quarks and gluons emerge, 

pdxfa{X,p) = 0, (3.5) 

with = {Ep,p). On integrating the constraint equation, Eq. ( |3.4| ), one has 
{Ep —jp — M'^)fa{X,p) = 0. The equation is the expression of the fact that the par- 
tons have to be on mass-shell, and is consistent with the quasiparticle assumption, 
Eqs. (|2.15|) - (|2.18| ) made in the first place. 

B. Fock self energies 

The next type of graph contributing to the mean field expansion is the Fock term. 
The generic diagrams for the quark and gluon self energies are shown in Fig. 30. 

Since we are particularly concerned with the collision integral, we start with the 
quark sector and examine as an example, the gain term generated by the Fock term 
I^p^{X,p). By inspection, one has 

^^-F:,iX,p) = I ^,S-^X,p,)G^-{X,p2){2nr6^'\p-p,+p,l 

(3.6) 



^Fig.3(b) is strictly speaking a vacuum polarization graph for the gluons 
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and Fpq is the Fock color factor for a single SU(Ai'c) group, 



The contribution to the collision term that this makes, using Eqs.( |2.19 ) and ( |2.20 ) 
is 

^Ccou = /^^ ^/S, + (X,p)(-^)5(po - E,)UX,pl (3.8) 

which, on inserting the explicit expressions for S ^(X, and {^^P) from 
Eqs.( 2.15|) and (|2.16|) leads to four distinct terms. 



TT TT 
El E2 



4 

x^T^E^^, (3-9) 



where 

Ti = 5{Ei -p',)SiE2-pl)fgiX,pi)f,iX,p2)f,{X,p) 

T2 = 5{Ei ~p\)5{E2+pl)UX,pi)UX, -p2)f,{X,p) 

T3 = +p°)/,(X, -p2)/,(X,p) 

T4 = -p^)/,(X,-pi)/3(X,p2)/,(X,p). (3.10) 

By attributing unbarred functions / to incoming particles and barred functions / to 
outgoing ones, one can see that T1..T4 correspond to the processes g — > qq, 0— > qqg, 
q —>■ qg and qg —>■ q. The last three of these are kinematically forbidden, while the 
former is possible, since the gluons are endowed with a finite mass. Performing the 
integrals over p^ and P2, Eq. (p.lO|) becomes 



yAM\l^j,(x,,,,:}Ux,rH)Ux,p). (3.11) 

The loss term is obtained in a similar fashion or by exchanging / with /, since the 
matrix element is symmetric. Combining both terms, the revised transport equation 
is obtained on integrating Eq.( |2.11| ) over the energy interval A+ as 
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J (27r)32^i (27r)32E2^ ^ 

(3.12) 

This is the final expression for the Fock transport equation. One can ahernatively 
introduce a cosmetic recombination or decay rate in which case Eq.( |3.12| ) can be 
written symbohcally as 



F 

qq^g 



x[f,iX,p2)f,iX,p,)f,{X,p) - /,(X,p2)/,(X,pi)/,(X,p)], (3.13) 
where F is the fiux factor, and 

..^^/.g^ (3.14, 

with the invariant phase space factor dQ given as dQ = (2it)^5^'^\p + pi — 

P2)d^P2/ii27Tf2E2). 

An analysis of the self energy graph 3(a)of the gluon sector, (X, p) along 

the previous lines leads to processes g — > gg, 0— ggg, and gg g, all of which are 
kinematically prohibited. One thus obtains 

,gain/loss)gluonic graph 

0. (3.15) 

This can be attributed to the fact that that the self energies are evaluated on shell, 
i.e. we may write 

^'pMaM^Po = Ep,p} = S^^-(,)(X,po = E„p) = 0, (3.16) 

which is the statement that an on shell particle cannot decay into two on shell 
particles of the same kind. 

The second graph in the gluonic case does not vanish. This self energy ^F~^(^b) that 
enters into the description of the gain in gluons, is precisely that given in Eq.( p.6| ), 
but with G^~{X,p2) replaced by ^^^(X, ^2)- The color factor in this case is also 
modified, being Fp^g = 1/25"". An analysis of the self energy along the same lines 
leads to the processes q —>■ qg, q qg, 0^ gqq and qq g, the last of which is the 
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only term that can contribute. Thus the time evolution of the gluon distribution 
function is given by 

x\M\l^^g[f,{X,p,)f,{X,p,)fg{X,p) - f,{X,P2)fg{X,p,)fg{X,p)]. 

(3.17) 

We conclude this section by commenting the result that while the Fock term 3(a) 
for gluons vanishes identically, the Fock term for the quark self energy does not. A 
term of this kind occurs in this model because the quarks are massless, while the 
gluons are massive. The relevance of this Fock term thus depends on the form of 
the underlying theory. In evaluating higher order graphs in the following section, 
it is useful to invoke the so-called "on-shell argument" , ignoring all terms of higher 
order of this kind that contain on-shell decays. This is necessary in order to render 
the problem tractable, but it is not strictly true, as has been indicated here. 

IV. THE COLLISION INTEGRAL - BEYOND THE MEAN FIELD 
A. Exact perturbative calculations that contain two loops 

Generic diagrams that contain two loops and which contribute to the quark and 
gluon sectors are shown in Fig. 4. In the quark sector, we denote these graphs as 
rainbow (R), ladder (L), cloud (C), exchange (E) and quark loop graphs (q-Loop), 
according to their topology. These graphs comprise all diagrams with two exchanged 
gluons. In the gluonic sector, there are more diagrams, as can be seen in the figure. 
Note that in this sector, a ladder-like diagram is topologically equivalent to the 
rainbow kind, and is therefore not included separately. Since we are specifically 
interested in the collisional term, we will concern ourselves only with the off-diagonal 
contributions that arise from these diagrams, as is required in Eq.( |2.11| ). In addition, 
these self energies are required on shell. 

To be specific, let us deal first with the quark sector, and consider the rainbow 
diagram first. Then there are four possibilities of ordering the lines in this graph, as 
is indicated in Fig. 5. However, invoking the on-shell argument, as explained in the 
previous section for the Fock diagram, only the first graph of this series can survive 
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while the other three vanish. Similar arguments can be applied to the remaining 
graphs, the L, C, E, and q-Loop diagrams. One can set up all possibilities of ordering 
the lines for each type of graph in Fig. 4, and determine which survive. The final 
set of all non-vanishing contributions can then be ascertained. These are shown in 
Fig. 6. 

We next note that the diagram of Fig. 6.4(b) is the exchange graph of Fig. 6.5, 
because both of these diagrams contain three off-diagonal quark Green functions. 
The similarity becomes more evident if one redraws the diagram of Fig. 6.4(b) as 
shown in Fig. 7. 

The self energies for two gluon exchange contribute to the quark collision term 

via 

4S=/ rfpoS(2)"+(X,p)5+-(X,p)- / rfpoS(2)+-(X,p)5-+(X,p) 

JA+ JA+ 

= -z^S(2)-+(X,po = i?p,p1/,(X,p-) +^^S(^)+-(X,po = Ep,p)f,{X,p), 
J^p tjp 

(4.1) 

where the off-diagonal self energies are constructed from the diagrams of Fig. 4 as 
where 

SSatk-quark(^, V) = ^^^^^ + ^2to"op(^, P) (4-3) 

and 

+ sgi;-(X,p) + sg)+-(X,p). (4.4) 

Here we have split up the combinations of the second exchange plus quark loop 
diagrams in Squark-quark, leaving the rainbow, cloud, ladder and first exchange graph 
to Squark-giuon- This subdivisiou in Eqs.( |4.3| ) and ( [4.4|) to J^Hi in Eq.( [4.1|) will be 
handled separately, as the first term will be seen to lead to elastic quark-quark and 
quark-antiquark differential scattering cross sections in the transport equation, while 
the Squark-giuon term will be seen to lead to processes involving gluons in the final 
state, such as the processes qq gg and qg ^ qg- 
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1. Quark-quark and quark- antiquark scattering cross sections. 



Explicit expressions for the quark loop and its exchange diagram self energies 
required in Eq.(|4.3|) are obtained as 



^(2)+- ^ 4 4E.2 /■ "Pi aP3 0^4 ,^ x4r(4)/ N 

x(27r)45W(p2 -P3 + P4) S+-(X,pi)G'++(X,p2) 

X (X, p,) (X, p4) G~ (X, P2) (4.5) 

and 

S,,, (X,p) = m ^i^y ^(^(^^(2-) \P-P.-P2) 

Xi27iy6^'\p2-P3+P4) S+-(X,pi)G++(X,p2) 

X 5+- (X, Pa) (X, p,)G-{X,p-ps), (4.6) 

where -Fq_Loop and Fe are color factors, that will be given explicitly in the following 
subsection. Since they do not affect our argument, we suppress them in the following. 

The integrated collision term, for loss, for example, from Eq.( [4.1|) can be directly 
evaluated, to give the quark loop and exchange contributions 



x{g++(X,p2)G-(X,p2) + G++(X,p2)G-(X,p-p3)} 

x(-^f )(-^f )(-^f )E^^> (4-7) 

-C/l -C/3 -C/4 j=i 



where 



Ti = 5(i?i + p?) 5iE, + pI) 5{E, + pI) f,{-p,) f,{-p,) fq{-p,) f^ip) 

T, = 6{E, + p1) 6{Es + pI) 6{E, - pi) f,{-p,) f,{-p,) f,{p,) f,{p} 

n = 5{E, + pI) 6{E, - pI) 6{E, + pI) f,{-p,) Up,) U-p^) f,{p) 

T, = 6{E, + p?) 5(^3 - pI) 6{E, - pI) /,(-pi) Up,) Upa) Up) 

n = S{E, - p?) 5iE, + pI) 5{E, + pI) Upi) U'Ps) U-Pa) Up) 

n = 5{E, - p\) 5{E, + pI) 5{E, - pI) Upi) U'Ps) fM Up) 

Tj = 6{E, - p\) 6{E, - pI) 6{E, + pi) Upi) IM) U'P^ Up) 

Ts = 5{E, - rf) 5iE, - p\) 5{E, - pi) Upi) Ups) Up^) Up)- (4i 
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One sees that there are eight terms, or eight processes in this expression. However, 
due to energy- momentum-conservation Ti,T2,T4,Tq and T7 vanish, leaving only T3, 
T5 and Tg. This is a direct consequence of the on-shell nature of the quasiparticle 
approximation. If this were relaxed, all terms would necessarily have to be included. 

We can reorganize this expression into a recognizable physical form by making 
some simple manipulations. Letting pi —pi for the antiquark states and perform- 
ing the p1, P3, p1 and the p2 integration by absorbing the appropriate 5-functions, 
we obtain 



TT 



y(2)loss ^ _ 4 4^ / ^-t^L ^' i^.', ^t"i /g y 

^ Ep7 (27r)32Ei (27r)'^2E3(27r)32E4^ ' 

X [8^^\p +P1-P3- Pa) fq{Pl) fq{P3) fq{P4) fq{p) 

x[G++{X,p + p^)G-{X,p + p,) + G++{X,p + p,)G-{X,p-ps)] 

+5^"^^ (P-Pl+P3- P4) fq(Pl) fqiPs) MPa) fq(p) 

x[G++iX,p-p,)G'-{X,p-p,) + G++iX,p-p,)G-{X,p + p3)] 

+ S^*^ {P-Pl-P3+ P4) fq{Pl) fqiPs) fqiPi) fq{p) 

X [G++ {X, p - p^)G- (X,p-p,) + G++ {X,p-p^)G-{X,p-p3)]}. 

(4.9) 

The first two terms of this expression can be combined if one makes the substitution 
Pi ^ Ps in the second term. The third term has a symmetry in pi and p^ and can 
be rewritten as one half the sum of two terms with pi and ps interchanged. The loss 
term then becomes 

.(2)loss _4.__4 ^ / d% 



•^coU.q y J (27r)32Ei (27r)32E3 (27r)32E4^ ' 



X [6^^^{p +P1-P3- Pa) fq{Pl) fqiPs) fq{PA) fq{p) 

X [G++ {X, p + p^)G- {X,p + p^) + G++ {X, p + p,)G-{X,p- ps) 
+G'++(X, p - ps)G- (X,p-ps) + p - ps)G- (X,p + p,)] 

+5^"^^ (P-Pl-P3+ Pa) fq{Pl) fqiPs) fq(PA) fq{p) 

X ^ (X, p - p,)G- {X,p-p,) + G++ (X, p - p,)G- (X, p - p3)] 

+G++(X,p -p3)G~(^,P -P3) + G++(X,p -p3)G~(^,P - Pi)]} . 

(4.10) 

Using the fact that [iG (p)]'^ = iG++(p) and making the substitution pi p4 in 
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the second term, one is able to identify the absolute values squared of the Green 

; has 



functions occuring in Jclii^°^^- One has 



J. 



(2)loss 



IT 



coll,q 



g m — 



Ep J (27r)32Ei (27r)32E3 (27r)32E4 



(27r)^(5(^)(p + Pi-p3-P4) 



+ 



- p,) + zG-{x,p- p,) Up) f,{p,) Up,) Up,) 
tG-{x,p + p,) + tG-{x,p- ps) ' Up) Up,) Un) Up, 



(4.11) 



Now one may recognize the scattering amplitude for elastic quark-quark scattering, 



iM 



qq^qq 



{pi 34) = {~igm)^[iG {p - p,) + iG {p - p,)] , 



and the scattering amplitude for quark-antiquark scattering. 



iM 



qq^qq 



{pl 34) = {-igm)''[iG~- (p + pi) + iG {p - p,)] , 



(4.12) 



(4.13) 



occurrmg m Eq. fPD , which may be concisely written as to give the final result 



7r 



t(2)1oss _ •■ 



d^i d% 



d^4 



Ep J {2tt)^2Ei (27r)32E3 (27r)32E4 



(27r)^5(^)(p + pi-j93-p4) 



X I2 l^^'^^iM ^ m Up) Upi) Ups) Up4) 
+ \M,^^,,ipi ^ 34)1' Up) Upi) Un) Up^)} 



(4.14) 



The Feynman graphs corresponding to these processes are shown in Figs. 8 and 
9 respectively. 



2. Quark-gluon scattering cross sections. 

We now turn our attention to the graphs of Fig. 6.1 to Fig. 6.4(a), which will 
lead to scattering processes that involve gluonic degrees of freedom. As in the 
previous section, the Feynman rules for non-equilibrium processes can be applied 
to these diagrams and the result Wigner transformed. This results in the following 
expressions for the self energies, 

x(27r)^5(^)(p2 -P3 -P4)G+-(X,pi)S++(X,p-p3) 
xG+-(X,p3) S+-{X,p,) S-{X,p- ps) (4.15) 
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for the rainbow diagram, 



xG+-{X,p,) S+-{X,p,) G-{X,p2) (4.16) 



for the ladder graph, 



x(27r)V^)(p2 -P3 -P4)G+-(X,pi)G±±(X,p2) 
xG+-(X,p3) 5+-(X,p4) ^^^(X,p-p3) (4.17) 

for the two cloud diagrams, and 

^em(^^P)- ^'^^^y (2vr)4(2vr)M2vr)M27r)4^^''^^ 

X (271)^5^(^2 - P3 - Pa) G'+-(X, Pi) 5++(X,p - ^4) 
xG+-(X,p3) 5+-(X,p4) ^-~(X,p - Pa) (4.18) 

for the first exchange diagram. Fr, F^, Fc and are appropriate color factors, that 
will be discussed in detail in the following subsection, but which will be suppressed 
here. Note that a factor 1/2 occurs in the expression for the ladder diagram because 
of the gluon loop. The expressions for are obtained from the ones for S*^^^^ 

by exchanging — and +. To construct the integrated collision term, one has as in 
the previous section to integrate over over an interval A±. which contains ±Ep{X) 
and one finds, using Eq. (|4.1| ) that the loss term of the collision incorporating the 
rainbow, cloud, ladder and first exchange graphs, is given as 

JcoL =ig m — ... ... ... ,^_,4 (27r) S"- >{p-pi- p2)5'- >{p2 - Ps - Pa) 



EpJ (27r)4(27r)4(27r)4(27r)4~ 
X (X, p - p3)S- (X, p-ps) + (X, p2)G- (X, P2) 

+G++{X,p2)S-iX,p- Pa) + G-(X,p2)S++(X,p - ps) 

+S++(X,p-p4)S-(X,p-p3)}(-.^)(-.^)(-.^)E^., (4.19) 

J J^l J^4 ■ 



^4 i=l 



where 
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S{E, + rf) 5(^3 + P^) 5{E, + pO) /.-(-pi) /g(-p3) /5(-P4) m 
S{E, + p\) 6{E, + pI) 6{E, - pI) f,{-p,) f,{-p,) Up,) Up) 
5{E, +pI) 5{Es-pI) S{E,+pI) f,{-p,) JM fsi-p,) Uip) 

5{E, + p?) 5{E^ - pI) 5{e, - pI) f,(-p,) Up,) Up) 
5{E, - p\) 6{E, + pI) S{E, + p°) Upi) U-ps) U-P4) Up) 
S{E, - p?) 6{Es + p°) 6{E, - pO) Upi) U-p^) IM) Up) 

S{E, - p?) S{Es - pI) S{E, + pI) Upi) Ups) U'P^) Up) 

5{E, - p?) 5{E^ - pI) 5{e, - pI) Upi) Ups) IM Up)- 



(4.20) 



Once again, eight terms result from this multiphcation. Now, again due to energy- 
momentum-conservation, Ti,T2,T3,T5 and Tg vanish, and we are left with three 
non- vanishing terms, T4, Tg and T7. 

Applying the same procedure as for Jclii"^^ as in the previous subsection, one 
can regroup the remaining terms to read 



X 



(2)loss 



coll,g 



4 4 

g m 



71 



I — 
J (277 



cPpA 



Ep J {27r)^2Ei (27r)32£;3 {27r)^2E, 



(27r)^5(^)(P + Pi-P3-P4) 



1 2 
- iG—{X,p + pi) + iS—{X,p - Pa) + iS—{X,p- Pi) 

^Up) UPi)UPs) fgiP4) 



+ 



iS— (X, p + Pi) + iG— {X, p - ps) + iS— {X, p - P4) 



xUp)Upi)Up3)Up4)}- 



(4.21) 



In order to identify the physical processes that give rise to these terms, we examine 
first all possible contributions to the annihilation process qq gg. The Feynman 
graphs for this within this model are shown in Fig. 10. The scattering amplitude 
associated therewith is 

- iMqq-,gg{pl 34) = {-igm)'^[iG~~ {p + pi) + iS~~ {p - P3) + iS~~ {p - p,)]. 

(4.22) 

In a similar manner, the elastic scattering process qg — > qg, which is shown in 
Fig. 11, has the scattering amplitude 

-iMqg^qgipl 34) = {-igm)^[iS~-{p + pi)+iG~-{p-ps)+iS—{p-p,)]. 

(4.23) 
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One can identify the absolute value squared of Eqs.( |4.22|) and (H:.23|) in Eq.( |4.21|) 

(2)loss 



and therefore J^Ji can be written as 



COlljg 



Ep J (27r)32Ei (27r)32E3 {2tt)^2E4 



X [I \M,,^,,{pl 34)1^ f,{p) f,{pr) Un) Up,) 

+ |-M„,(pi ^ 34)1^ f,{p) Up,) Up,) Up,)] . 

The complete loss term is obtained by adding Eq. (^4. 14|) and ( [4 .241) , 

^(2)loss _ ^(2)loss y-(2)loss 
"^coU ~ '^coU.q + ^coll.g ■ 



(4.24) 



(4.25) 



The gain term can be constructed by replacing / ^ / in the complete loss term. 
With the relation 

|XP dQ 



da 



dn \vp - Vi\2Ep2Ei dn 



(4.26) 



and the phase space factor 



g = (27r)V4)(p + Pi-p3-P4)- 



d^p, d^p. 



(4.27) 



(27r)32E3 {27t)^2E,' 

the final form for the Boltzmann equation, calculated to second order in the number 
of exchanged gluons, is for quarks {a = q) 

d^Pi 



2pdxUX,p) = J - vi\2Ep2E, 



daj 



i=i 



ab^cd 



(27r)32Ei 

UP)fb{Pl)fc{P3)fd{p4) - fa{p)fb{pi)fc{n)fd{Pi] 



(4.28) 



where partons 6, c, and d can be a quark, antiquark or gluon, and j labels the four 
processes j = 1...4 corresponding to qq —>■ gg, qg qg , qq qq and qq qq. The 
Sj are symmetry factors si = S3 = 1/2 and S2 = = 1. 

The transport equation for gluons can be obtained in an analogous way and 
calculated for two loops, it takes the same form as Eq.( |4.2(j| ) with a = g. Then j 
labels the four processes j = 1...4 corresponding to gg gg, gg qq, gq — > gq and 
gq — > gq. The appropriate symmetry factors are Si = 1/2 and S2 = S3 = S4 = 1. 
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B. Color expansion 



We now deal with the color factors, which we have neglected so far. We calculate 
them for one SU (N) color group. The overall color factor for both color groups is 
then obtained by squaring it. 

The matrices {t°')ij are the matrices of the color group in the representation of 
the quarks, while {T"-)bc = —ifabc are the color matrices in the adjoint representation 
and fabc are the structure constants of the color group, see Eq.(|2.2|). The t"'s are 
normalized to 

tiiet') = ^5ab. (4.29) 

The "square" of the generator in some representation must be proportional to the 
unit operator (Schur's Lemma). Therefore 

{ndnjk = CF6,k (4.30) 

and 

TbdTdc = fbadfcad = Ca^c-, (4-31) 

where the numbers Cp and Ca are the Casimir operators of the fundamental and 
adjoint representation, respectively. They take the values (see for example pT| ) 

N'^ - 1 

Cf = (4.32) 



2N, 



c 



and 

Ca = N,. (4.33) 

Consider now the quark self energies that were evaluated in the previous subsection. 
Let i denote the external parton index. It is therefore not to be summed over. The 
color factor for the rainbow graph is 

b a a b 2 (^c - 1)^ 

Fr = tijf^kttitii = Cp6ii = — — 5ii. (4.34) 
For the ladder graph, one finds 

-1 

Fl = i-ifabc){-ifcbd)tijtji = CA^adttjtji = CaCf^U = Sii- (4.35) 
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For the cloud graph, one obtains 



Fc = {-ifacb)t%t%tli = - "\ ^ Sii, (4.36) 



where the relation (see for example ||2l[| ) 



tfabcf't' = (4.37) 



has been used. The color factor for the exchange graph is 



a b a b - 1 

Fe = ttjtjA'tu = 4/Y2~^"' (4.38) 



where the relation 



t^th^ = (4.39) 



has been used. Finally, for the quark loop graph the color factor is given by 

- 1 

Fq-Loop = tt/.Mf't') = -J^S^^- (4-40) 

Up to this point, we have made a semiclassical expansion that involves keeping only 
the leading term in expanding the exponential in Eg. ( p.l4| ) (here the factor h has 
been set to one.) In addition, we have examined sets of diagrams organized according 
to the number of interaction lines, i.e. according to the coupling strength. We have 
found that all graphs in a class are required in order to build up the cross sections 
that ultimately occur in a Boltzmann like equation. However, at the level of two 
exchanged gluons, we are already faced with five types of graphs, and this number 
increases rapidly with the number of exchanged gluons. One possible simplifying 
assumption is the additional imposition of an expansion in the inverse number of 
colors. According to such a criterion, the ladder, the rainbow and the cloud diagrams 
are leading, since their color factors for one color group are of order 0{N^) while 
for the quark- loop diagram it goes as Nc and for the exchange diagram only as A'"^. 

Since the ladder and the rainbow diagram lead to cross sections involving gluons 
while the quark-loop diagram leads to elastic quark- ( ant i) quark cross sections, one 
can conclude that the quark degrees of freedom are suppressed in comparison with 
the gluon degrees of freedom. This is in agreement with the results of an evaluation 



of the quark-quark scattering amplitude within this model |T^, in which the quark 
degrees of freedom are neglected, however due to kinematical reasons. 
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Although the ladder and the rainbow diagram are both of order 0{N^), the ratio 
of their color factors for one color group is not one, but 

which is 4/9 1/2 for Nc = 3. Since in the rainbow diagram the second gluon 
couples at the quark line while in the ladder diagram it couples at the first gluon, 
two quark-quark-gluon vertices are suppressed by a factor 4/9 per color group in 
comparison with two 3-gluon vertices. Thus, there is no strict ordering of the gluon 
graphs according to a single class of diagrams, in an expansion in l/N^. Although 
the ladder graphs and the processes that they lead to appear largest, one should 
not that the symmetry factors of the other graphs compensate for this. A numerical 
study is essential to determine the actual order of magnitude of each graph. 

As a final comment, note that an expansion in color also incorporates the coupling 
strength. Assuming that g ~ 1/iVc) we find that the Fock term ~ g^N^ and the 
ladder diagram ~ g'^N^ are of the same order. 



C. Three parton production 

Since many diagrams contribute to the self energy with three exchanged gluons, 
we again make an additional expansion in the inverse number of colors as explained 
in the previous Subsection, and consider only diagrams that do not contain quark 
loops. In the case of two exchanged gluons, we have seen that the ladder and 
the rainbow diagrams give rise to the squared scattering amplitudes of individual 
channels. Therefore, we conclude (and will check later) that for the case with three 
exchanged gluons, the three-rung ladder diagram, the rainbow diagram with three 
bows and all possible mixtures between them will lead to the squared scattering 
amplitudes of individual channels of processes involving three gluons. The generic 
form of these diagrams is given in Fig. 12. (Note: diagrams that would generate 
crossed terms between the single channels are shown in Appendix C). 

A direct evaluation of the self energy diagrams in Fig. 12 should give rise to 
cross sections of processes that involve three gluons, i.e. one expects, for example, 
the three gluon production qq ggg but also multiparton production such as 
qg qgg. The generic Feynman diagrams for these processes are shown in Figs. 13 
and 14 respectively. 
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To be specific, let us consider now the ladder diagram of Fig. 12(a), in particular 
Tr^ [X^p). There are sixteen possibilities of arranging the signs that can be 
associated with the vertices of the first and second rungs. Invoking the on-shell 
argument and energy-momentum conservation, the possibilities become restricted 
to the single diagram shown in Fig. 15. The explicit expression for this ladder 
diagram is 

(3)+_ 6 6 [ dj)! d% d% d% d% d% 

^ ^ ~ '2^ 7 (27r)4 (27r)4 (27r)4 (27r)4 (27r)4 {2it)* 

x(27r)V^)(p -pi -P2)(27r)^5(^)(p2 -P3 -P4)(2vr)V^)(p4 -P5 -Pe) 

X S+- (X, p^) G++ (X, p2)G+- (X, ^3) G++ (X, p,) 

X G+- (X, p,) G+- (X, pe)G" (X, p,)G- (X, P2). (4.42) 

Once again, one may insert the quasiparticle ansdtze for the quark and gluonic 
Green functions from Eqs.( |2.15|) - (|2.18|) , and construct the contribution of this 
ladder term to JcUl"^^ as in Subsection IV A by taking its on-shell value 

= -^^S(=^)+-(X,po = Ep,p)f,{X,p). (4.43) 

Multiplying out the product of Green functions in terms of the quark and gluonic 
distribution functions enables us to identify all processes that can occur. The details 
of this cumbersome calculation are given in Appendix B. Here we simply note that 
the constraints of energy-momentum conservation admit all processes which have at 
least two partons in both the initial and final state. That is to say, in addition to the 
parton production processes qq — > ggg and qg qgg, that were already mentioned, 
one also obtains the matrix elements for the processes qqg gg and qgg —* qg. 

In order to complete the calculation, the remaining diagrams of Fig. 12 and 
additional self energy graphs that lead to mixed terms of the scattering amplitudes 
must be calculated in an analogous way. There are very many such diagrams that 
build the mixed terms between the scattering amplitudes that were shown in Figs. 13 
and 14. For completeness, we show these in Appendix C. We will not however 
calculate these graphs in any detail here. 

The final form of the Boltzmann equation, calculated to third order in the number 
of exchanged gluons, reads 

o fi f fv ^ f ^^P^ ^^^2 d^p3 4 
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X < ip + pi-p2- Ps) -^i \-Mj{qa -> bc)\' 



X [fq{P)fa{Pl)fb{P:i)fc{Pi) - fq{P)fa{Pl)fb{n)fc{PA 

+ / ('f\loT? I ^^^^ iP + Pi-P2-P3- Pa) Y1 *i bcd)f 

X fq{P)fa{Pl)fb{n)fc{Pi)fd{n) - fq{P)fa{Pl)fh{p?,)fc{PA)Id{Pb) 

8 

+5^'^\p + P1+P2-P3- Pi) Sj \Mj{qab cd)f 

i=7 

X fq{P)fa{Pl)fh{P3)fc{Pi)fd{Pb) - fq{P)fa{Pl)fb{P3)fc{Pi)fd{P5)] } , 

(4.44) 

where a, b, c and d are partons and j = 1...4 labels the four processes given after 
Eq.( [4.28[ ) while j = 5. ..8 labels the four processes corresponding to qq — > ggg, 
(19 ~^ 19 9 1 QQ9 ~^ 99 ^-iid qgg qg. The symmetry factors for j = 5... 8 are given as 
S5 = 1/3! and Sg = = Sg = 1/2. 



D. n parton production 

Of the self energy diagrams of order 0{g^^), the ladder diagram, the rainbow 
diagram and all possible mixtures between these two are leading in an expansion in 
l/Nc- On evaluation, these diagrams lead to the scattering process qq ng and all 
possible crossed processes, such as qqg {n — l)g, qg q{n — l)g, ... in which at 
least two partons occur both in the initial and final states. The Boltzmann equation 
for quarks then reads 

d^k d^pi _d^pm+n 

P m,n=l ' 



pdxfq{X,p} = - E /- 



(27r)32Efc (27r)32Ei {27r)^2Em+n 

X [d^'^'ip + k+pi + ... +Pm-1 -Pm- ••• -Pm+n) 

xSm-iSn+1 \M{qq{m - l)g ^ (n + l)g)f 



(2vr) 



X 



fq{P)fq{k)fg{pi)--fg{Pm^l)fg{Pm)--fg{Pm+n 
-fq{P)fq{k)fg{Pl)--fg{Pm-l)fg{Pm)--fg{P m+n) 
+ 5^^\p + pi + ... +Pm-k -pm+1 - ••• -Pm+n)SmSn \M{qmg 
X fq{p)fg{Pl)--Jg{Pm)fq{k)fg{prn+l)---fg{Pm+n) 



qng)\' 
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-fg{p)fg{Pl)-f9{Pm)fq{k)fg{pm+l)-fg{Pm+n) 



(4.45) 



with the symmetry factors s„ = 1/?^!, which is our final result. 

V. SUMMARY AND CONCLUSIONS 

In this paper, we have worked out the transport theory for the scalar partonic 
model. We have placed particular emphasis on the collision term examining the 
collision integral first in the mean field (Hartree and Fock approximations) and 
then including more loops. On making an analysis of the color structure, we have 
found that a 1/Nc expansion leads to a suppression of the quark loops and their 
exchange diagrams, much as occurs in QCD. Thus it turns out that all remaining 
self energy diagrams that lead to gluon production and quark-gluon scattering are 
of the same order. We have not been able to find any convincing argument to 
favor an expansion that only includes ladder graphs, as is for example the dominant 
contribution in elastic quark-quark scattering in the Regge limit, when considered 
on its own within this model. 

This can simply be traced back to the fact that the uppermost line in the quark 
or gluon self energy is closed and therefore possesses a momentum that must be 
integrated over. This fact reemerges in the Boltzmann equation in that all momenta 
except for that of the particle under study, are integrated over. An additional 
assumption, e.g. that the uppermost gluonic momentum should take on a specific 
value, is required in order to generate a leading logarithmic expansion. 

We have also analysed the situation involving three loops, which leads to three 
parton production, and find that initial state interactions also enter into the pic- 
ture. Finally, a generalized Boltzmann like equation that includes multiple gluon 
production, but which also incorporates initial state interactions, is given. 

In concluding, we comment that further study is required in order to establish 
the factors that could lead to mass renormalization in this model and which has not 
yet been addressed. 
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APPENDIX A: GREEN FUNCTIONS, TRANSPORT AND CONSTRAINT 

EQUATIONS 

In this appendix, we give a brief guideline for the derivation of the transport and 
constraint equations, Eqs.(p.8|) and ( |2.9| ). 

The Schwinger-Keldysh Green functions defined in Eqs.( |2.5| ) and ( |2.6| ) can be 
summarized in a compact matrix form 

using the generic notation already introduced in Section JIB. From the definition of 
the Green functions, the relation 

D--{x, y) + y) = D~+{x, y) + D^'{x, y) (A2) 

follows, showing that the four components D^^ are not independent. We define the 
retarded and advanced Green functions in the standard way as 

:= Q{xo- y,)[D+-{x,y) - D-+{x,y)] 

= D-{x, y) - D-+{x, y) = D+'{x, y) - D++{x, y) (A3) 
D^{x,y) := -e(yo - xo)[D+-{x,y) - D-+{x,y)] 

= D-{x,y) - D+-{x,y) = D-+{x,y) - D++{x,y). (A4) 

The equations of motion that the Green functions satisfy are 

+ M^)D{x,y) = -a,6^'\x - v) + J d^z a,U.{x, z)D{z,y), (A5) 

given in terms of the irreducible proper self energy 
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and 



'1 0' 
,0 -1, 



(A7) 



In Eg. (|A5|) , M is the free bosonic parton mass. 

The four components of the self energy are also not independent. From their 
definition, the relation 



(A8) 



can be seen to hold. The retarded and advanced self energies are defined to be 

U''{x,y) = U-{x,y) + U-^{x,y) 

n^{x,y) = n-{x,y) + n+-ix,y). (A9) 
We now consider specifically the equation of motion for D ^. This reads 

+ M^)D-+{x, y)=J d'z{U-{x, z)D-+{z, y) + n-+(a;, y)} 

= J d'^z{U^{x,z)D-+{z,y) -U+-{x, z)D-+{z,y) 

+n-+(a;, z)D+-{z, y) - I[-+{x, z)D^{z, y)} (AlO) 

while the conjugate equation is 

+ M2)D-+(x, y) = -J d^'ziD-^ix, z)n++(;2, y) + D—{x, z)Il~+{z, y)} 

= J d'z{D-^{x,z)U^{z,y) - D^{x,z)U~^{z,y)}. (All) 

Now a Wigner transform of both equations is performed to yield 
■1 



-□ 



X 



X 



U'^{X,p)AD-+{X,p) - Il+-{X,p)AD-+{X,p) 
U~+{X,p)AD+~{X,p) - U'+{X,p)AD^{X,p) 



(A12) 



and 



-Bx + ipdx - p"^ + M"^ 



D-+(X,p) = 
5-+(X,p)AS^(X,p) - ^«(X,p)AS-+(X,p), 



(A13) 
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with the differential operator 



A:=exp\ — {^dxdp- dpdx 



Subtracting Eg. ( |A12|) from Eq. ([A13D gives the so-called transport equation, 



(A14) 



(A15) 



while adding them yields the so-called constraint equation, 

(^□x - 2p' + 2M2) D-+{X,p) = /+, (A16) 
that were quoted as Eqs.(^.8D and (^.91), and is as given in Eqs.( |2.1()D to ( |2.13| ). 



APPENDIX B: THREE PARTON PRODUCTION 

In this Appendix, the explicit evaluation of the ladder graph for three parton 
production is given. On inserting the quasiparticle Green functions from Eqs.( ^.l5| ) 
to (|2.18|) in Eq. ([4.42|) and after performing the p2, Pi, Pi, p^, P5 and pg integration 
by absorbing the appropriate 5-functions we obtain from Eq. (|4.43|) the result 



t(3)1oss 
^coU 



1 g g TT /• d^pi d^p3 d^p5 d^pe 

Q JYl / 

2^ Ep J (27r)32Ei (27r)32E3 {2-Kf2E^ {2-nf2EQ 



XG++{P3 +P5+ Pa)G -{P3 +P5+ P6)G'++(p5 + P6)G—{P5 + P&) 

16 

X {2Tx f5^^ (p _ _ p3 _ p5 _ pg) ^ Ti, (Bl) 



i=l 



where the Tj are given as 



Tl = fq{p)fq{Pl)fg{n)fg{Pb)fg{P&) 
T2 = fq{p)fq{Pl)fg{.n)fg{P5)fg{-n) 
T3 = fq{p)fq{Pl)fg{P3)fg{-P5)fg{P6) 
Ta = fq{p)fq{Pl)fg{n)fg{-Pb)fg{-P&) 
T5 = fq{p)fq{Pl)fg{-n)fg{Pb)fg{P&) 
Tq = fq{p)fq{Pl)fg{-P3)fg{P5)fg{-m) 
Tj = fq{p)fq{Pl)fg{-n)fg{-Ph)fg{P6) 
Ts = fq{p)fq{Pl)fg{-p-i)fg{-Ph)fg{-P&) 
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~* \ c / ~* \ c / ~* \ c / ~* \ 

-Pl)fg{P3)f9{P5)fg{P6) 




= fMfgi- 


-PijfgiPsjfgiPrdfgi-Pa) 


Tu 


= fMM- 


-Pl)fgm)fg{-Pb)fg{P(i) 




= f<,iP)fgi- 


-Pl)fg{P3)fg{-P5)fg{-P6) 




= mM- 


-Pl)fg{-n)fg{Pb)fg{PQ) 


Tu 




-Pl)fg{-P3)fg{P5)fg{-P6) 




= mui- 


-Pl)fg{-n)fg{-P5)Ig{P%) 


Tl6 




-Pi ) fg ( -P3)fg ( -P5 ) fg ( -Ps) 



By attributing unbarred functions / to incoming particles and barred functions 
/ to outgoing ones, one can make the identification of Tig corresponds with the 
annihilation process qqggg -^0, Ti with q —>■ qggg, Tg with qggg q while T12, Tu 
and Ti5 correspond to the process qqgg — > g. These processes are all kinematically 
forbidden. Only processes with at least two partons in both the initial and the final 
state are kinematically allowed: these are T2, T^ and T5, which lead to qg qgg, 
T4, Tq and T7, which lead to qgg — > qg, Tiq, Tu and T12, corresponding to qqg — > gg, 
and finally Tg, which gives rise to qq — > ggg. 



APPENDIX C: REMAINING SELF ENERGY DIAGRAMS 

In Section IIIC, we gave the Feynman graphs which lead to the square of the 
individual scattering amplitudes. There are many diagrams which give rise to the 
mixed terms that also are necessary when forming the square of the full amplitude. 
For completeness, we indicate which generic diagrams also contribute to the self 
energy E^^^ and lead to these mixed terms between the scattering amplitudes of the 
individual channels qq — > ggg, qg — > qgg, qqg — > gg and qgg — > qg. These are given 
in Fig. 16. The last six diagrams are not symmetric, and it is to be understood that 
the mirror reflected diagrams must also be taken into account. 



30 



REFERENCES 

[1] See for example, C.-Y. Wong, Introduction to High-Energy Heavy-Ion Collisions, 
(World Scientific, Singapore, 1994). 

[2] K. Geiger and B. Miiller, Nucl. Phys. B369, 600 (1992). 

[3] K. Geiger and B. Miiller, Phys. Rev. D 50, 337 (1994). 

[4] K. Geiger, Phys. Rev. D 54, 949 (1996). 

[5] U. Heinz, Phys. Rev. Lett. 51, 351 (1983); H.-T. Elze and U. Heinz, Phys. Rep. 
183, 81 (1989). 

[6] H.-T. Elzc, M. Gyulassy, and D. Vasak, Phys. Lett. B177 402 (1986); Nucl. 
Phys. 276, 706 (1986). 

[7] W. Botermans and R. Malfliet, Phys. Rep. 198, 115 (1990). 

[8] K.-C. Chou, Z.-B. Su, B.-L. Hao, and L. Yu, Phys. Rep. 118, 1 (1985). 

[9] L.D. Landau and E.M. Lifschitz, Physikalische Kinetik (Akademie Verlag, 
Berlin, 1986), Vol.10. 

[10] S. Mrowczyhski and U. Heinz, Ann. Phys. (N.Y.) 229 , 1 (1994). 

[11] J.E. Davis and R.J. Perry, Phys. Rev. C 43, 1893 (1991). 

[12] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 5 (1961); ibid 124, 246 (1961). 

[13] S.P. Klevansky, A. Ogura, and J. Hiifner, Ann. Phys. (N.Y.) 261, 261 (1997). 

[14] J.C. Polkinghorne, J. Math. Phys. 4, 503 (1963); 4, 1393 (1963); 4, 1396 (1963). 

[15] J.R. Forshaw and D.A. Ross, Quantum Chromodynamics and the Pomeron 
(Cambridge University Press, Cambridge, 1997). 

[16] U. Heinz, Phys. Rev. Lett. 51, 94 (1986). 

[17] P.A. Henning, Nucl. Phys. A546, 653 (1992). 

[18] P. Rehberg, Phys. Rev. C 57, 3299 (1998). 

[19] J. Schwinger, J. Math. Phys. 2 , 407 (1961); L.V. Keldysh, JETP 20, 1018 

31 



(1965). 

[20] J.D. Bjorken and S.D. Drell, Relativistic quantum mechanics (Mc Graw-Hill, 
New York, 1964). 

[21] Yu.L. Dokshitzer, in Lecture Notes in Physics, edited by F. Lenz, H. Griefiham- 
mer, and D. StoU ( Springer- Verlag, Heidelberg, 1997), Vol.496, p. 87. 



FIGURES 
Fig. 1. Closed time path. 

Fig. 2. Quark and gluon generic Hartree self energies. Solid lines refer to quarks, 
wavy lines to gluons. 

Fig. 3. Quark and gluon generic one loop self energies. The quark self energy plus 
the first gluon self energy are Fock diagrams, while (b) of this figure is a polarization 
insertion. 

Fig. 4. Generic diagrams for the quark and gluon self energies that contain two 
loops. 

Fig. 5. The rainbow graphs. 

Fig. 6. All diagrams contributing to E(^)+~. 

Fig. 7. Diagram of Fig. 6.4b. 

Fig. 8. The t and u channel Feynman graphs for elastic quark-quark scattering. 
Fig. 9. The s and t channel Feynman graphs for elastic quark- antiquark scattering. 
Fig. 10. The s, t and u channel Feynman graphs for the process qq — > gg. 
Fig. 11. The s, t and u channel Feynman graphs for the process qg ^ qg- 
Fig. 12. Generic diagrams for the quark self energy that contain three exchanged 
gluons and which give rise to squared scattering amplitudes for individual channels. 
Fig. 13. Feynman graphs for the process qq — > ggg. 
Fig. 14. Feynman graphs for the process qg — > qgg. 

Fig. 15. The only three rung ladder contribution to the quark self energy. 

Fig. 16. Generic graphs for the quark self energy that are required for constructing 

the mixed amplitudes that lead to scatterings with three gluons. 
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